High-energy pseudogap in degenerate Hubbard model induced via Hund coupling 
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Hund coupling in the degenerate five-band Hubbard model near n = 6 occupancy is shown to give 
rise to a significant depletion of spectral weight above the Fermi level. Calculations within dynamical 
mean field theory combined with exact diagonalization reveal that this pseudogap is associated with 
a collective mode in the self-energy caused by spin fluctuations. The pseudogap is remarkably stable 
over a wide range of Coulomb and exchange energies, but disappears for weak Hund coupling. The 
implications of this phenomenon for optical spectra of iron pnictides are discussed. 

PACS. 71. 20. Be Transition metals and alloys - 71.27+a Strongly correlated electron systems 



The discovery of superconductivity in iron pnictides 
[l[ has stimulated intense discussions concerning the role 
of correlation effects in these compounds. In contrast to 
high-Tf, cuprates, which have antiferromagnetic Mott in- 
sulators as parent compounds, pnictides are correlated 
magnetic metals that tend to show significant devia- 
tions from Fermi-liquid behavior. Moreover, as a result 
of the multi-band nature of pnictides, the interplay of 
Coulomb and exchange interactions give rise to phenom- 
ena not found in cuprates. The importance of Hund cou- 
pling in pnictides and chalcogenides was recently pointed 
out in several theoretical [2-l6| and experimental stud- 
ies Optical data on paramagnetic LaFePO 
and BaFe2As2 0-0] reveal a high-energy pseudogap not 
compatible with normal metal behavior. This pseudo- 
gap differs fundamentally from the low-energy gap in the 
antiferromagnetic spin-density wave phase. Also, angle- 
resolved photoemission spectra for doped BaFe2As2 ex- 
hibit a depletion of spectral weight near the Fermi level 
that differs from the superconducting gap 

To investigate the influence of Coulomb correlations on 
the electronic properties of iron pnictides, various groups 
[1, @, ll^42l| have used dynamical mean field theory 22 1 
(DMFT) combined with single-particle Hamiltonians de- 
rived from density-functional theory. These calculations 
typically revealed moderate to strong effective mass en- 
hancement, in agreement with experimental data. Corre- 
lations were also shown to give rise to nonzero low-energy 
scattering rates [l^, 13 [3], indicating bad metallicity. 
On the other hand, as a result of the complex band struc- 
ture of pnictides, all five d bands are important, so that 
the many-body features exhibit a marked orbital depen- 
dence. Moreover, the Fe 3d density of states generally 
shows two main features separated by a deep minimum 
slightly above the Fermi energy. In view of these multi- 
band characteristics, it is difficult to distinguish genuine 
many-body features from single-particle properties. For 
instance, it is presently not clear to what extent pseudo- 
gaps in the interacting density of states are induced by 
Coulomb correlations or band structure effects. 

The aim of the present work is to unravel these com- 
peting influences. For this purpose we have performed 



DMFT calculations for a simplified Hamiltonian consist- 
ing of five degenerate semi-elliptical bands. Since the 
density of states is featureless, correlation induced spec- 
tral modifications are easily identified. In addition, by 
scanning a wide range of Coulomb and exchange ener- 
gies, it is feasible to investigate their respective roles. 
The main result of this work is that many-body effects are 
much more sensitive to the magnitude of Hund's rule cou- 
pling J than to the intra-orbital Coulomb repulsion U . In 
particular, exchange interactions give rise to a pseudogap 
above Ep, which persists in the full range of realistic val- 
ues of J, but disappears at small J. Within the present 
model, this pseudogap can be linked to a resonance in the 
self-energy caused by spin fluctuations. We also argue 
that these correlation effects can be understood qualita- 
tively by viewing the system at n = 6 occupancy as a 
doped n = 5 Mott insulator. Because of the high orbital 
degeneracy, the paramagnetic half-filled system is metal- 
lic at small J, but insulating already at moderate Hund 
coupling. Upon electron doping, a narrow quasi-particle 
peak appears in the density of states, with a maximum 
below and a pseudogap above Ep. Thus, as a function 
of J, the pseudogap at n = 6 exists only in proximity to 
the Mott phase of the half-filled system. 

To motivate the present work, we show first in Fig. 1 
a comparison of the electron spectral distribution of 
FeAsLaO derived within DMFT with the bare density 
of states obtained from the effective id single-particle 
Hamiltonian provided in Ref. |23|. Both spectra exhibit a 
striking minimum slightly above Ep. As impurity solver 
we use exact diagonalization (ED) [i^l- Each d orbital 
hybridizes with one or two bath levels, giving total clus- 
ter sizes ris = 10 or Us = 15, respectively. The interact- 
ing Hamiltonian and additional details concerning the 
ED DMFT approach are specified in Ref. Coulomb 
and exchange interactions are treated in a rotationally in- 
variant manner. The temperature of the Matsubara grid 
corresponds to T « 0.01 eV, while the interacting cluster 
Green's function is evaluated at T = for computational 
reasons. The spectral distribution at real ui is obtained 
by using the routine ratint to extrapolate the lattice 
Green's function from the imaginary axis. Typically, a 
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FIG. 1: (Color online) Fe 3d quasiparticle spectra of 
FeAsLaO, calculated within ED DMFT (n« = 15) for U = 
3 eV, J = 0.75 eV (solid red curve). The bare 3d density of 
states is indicated by the dashed blue curve :23] . 



few hundred Matsubara points are used in this extrap- 
olation, with a frequency dependent broadening. This 
procedure yields reliable results at low frequencies. Un- 
certainties arise mainly at higher frequencies and affect 
the position and width of the Hubbard bands. 

The quasiparticle spectrum of FeAsLaO is seen to ex- 
hibit two main peaks, below and above Ep, separated by 
a narrow pseudogap slightly above Ep, and a lower Hub- 
bard band below the Fe 3c? band region. The bare den- 
sity of states also exhibits pronounced peaks below and 
above Ep, and a deep minimum about 0.3 eV above the 
Fermi energy. Comparing these spectra one is tempted 
to conclude that Coulomb interactions lead to a signifi- 
cant narrowing of the occupied 3d bands and to a strong 
shift towards Ep. Only a small fraction of the occupied 
spectral weight 0.15/0.6 per spin band) remains in the 
coherent peak near Ep, while the main part 0.45/0.6) 
is distributed over the wide incoherent region. In the un- 
occupied region, the effect of Coulomb interactions is less 
severe. The main density of states peak is merely broad- 
ened and the minimum about 0.3 eV above Ep is approx- 
imately preserved. Similar quasi-particl e sp ectra have 
been obtained for a variety of pnictides (6l,ll3l- 21. 26. 27|. 
In principle, integer occupancy n = 6 can also sustain a 



Mott phase. This phase occurs, however, at much larger 
Coulomb energies [4, 16, Is! ]. 

Although the scenario discussed above appears plau- 
sible, it is nontrivial to disentangle genuine many-body 
features from the complex single-particle aspects of the 
density of states. In order to identify interaction effects 
due to U and J we have carried out ED DMFT calcula- 
tions for a Hamiltonian consisting of five identical sub- 
bands with semi-elliptical density of states, where the 
width W ~ 4 cV corresponds to typical pnictide com- 
pounds. Fig. 2 shows spectral distributions for a wide 
range of Coulomb and exchange energies, [/ = 3, . . . , 5 eV 
and J = 0.5, . . . , 1.0 eV. Most previous DMFT studies of 
pnictides used U and J values within this range. For 
simplicity, the present calculations are performed using 
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FIG. 2: (Color online) Low-energy region of spectral distri- 
bution of degenerate five-band model at occupancy n = 6, 
calculated within ED DMFT {us — 10) for several Coulomb 
energies: U — 3 eV (solid red curve), f/ = 4 eV (short-dashed 
blue curve), and [7 = 5 eV (long-dashed green curve), (a) 
J = 0.5 eV, (b) J — 1.0 eV. Spectra at intermediate values 
of J are similar. The upper panel also shows the results for 
U — 3 eV, J = 0. The bare density of states is indicated by 
the dotted curve. 



one bath level per 3d orbital {ug = 10). The results agree 
qualitatively with more accurate ones for two bath levels 
per orbital. Even though U and J span a wide range, the 
spectra are seen to be remarkably similar. In all cases, 
the interacting density of states exhibits peaks below and 
above Ep, separated by a pseudogap slightly above Ep. 
The peak-dip structure near Ep is reminiscent of the 
one in the quasi-particle spectrum shown in Fig. 1. Since 
in the present model the bare density of states is fea- 
tureless, the pseudogap is entirely due to the frequency- 
dependent self-energy. Moreover, the pseudogap depends 
only weakly on the values of U and J within the range 
quoted above, but it disappears at small J, so that only a 
broad quasiparticle peak remains. An example is shown 
in Fig. 2(a) for [/ = 3 eV, J = 0. These results suggest 
that the pseudogap above Ep is a generic feature caused 
by multiband Coulomb correlations within the 3d shell 
and that its existence depends crucially on realistic val- 
ues of Hund coupling. The paramagnetic quasi-particle 
distribution of actual pnictides should therefore consist 
of a combination of correlation features associated with 
J and signatures related to the bare density of states. 

To illustrate the effect of Coulomb correlations on the 
self-energy, we show in Fig. 3 the imaginary part of 
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FIG. 3: (Color online) (a) Imaginary part of self-energy for 
degenerate five-band model, calculated within ED DMFT 
(ris = 15) as a function of Matsubara frequency for vari- 
ous occupancies ai U — 4 eV, J = 0.75 eV. (b) Variation 
of self-energy (left) and spin susceptibility (right) with Hund 
coupling for [/ = 4 eV, n = 6 (ris = 10) . 

E(jcj„) for several occupancies. For n > 6.2, Im I](ia;„) ^ 
cjn at low energies, so that the system is a correlated 
Fermi liquid with a doping dependent effective mass en- 
hancement m*/m w 3,..., 6. At smaller occupancy, 
ImE reveals a nonzero onset, leading to bad- metallic be- 
havior. As shown previously [l9| . the spin correlation 
function then changes from Pauli to Curie Weiss behav- 
ior, indicating a spin freezing transition, in close analogy 
to the one found in the degenerate three-band model near 
n = 2 [2^. In addition, ImE(ia;„) exhibits a novel kink 
near cij„ « 0.1, . . . ,0.2 eV, which is weak in the Fermi- 
liquid region, but becomes more intense at lower occu- 
pancy. As shown in panel (b) , this kink is related to a res- 
onance in ImE(aj), which is the origin of the pseudogap 
above Ep in the spectra shown in Fig. 2. Remarkably, the 
resonance disappears at small J. For J = 0.5, . . . , 1.0 eV, 
the local spin correlation function reveals a maximum at 
about 0.1 eV, suggesting that the resonance in ImE(LLi) 
corresponds to a collective mode induced by spin fluc- 
tuations associated with Hund coupling. Thus, states in 
the resonance region have a greatly reduced lifetime. Ac- 
cording to Kramers-Kronig relations, ReE(a;) exhibits a 
positive slope near the resonance, so that spectral weight 
is removed from the pseudogap region. Beyond this re- 
gion, the slope of ReE(a;) becomes again negative, giving 
rise to a kink in the dispersion of energy bands |29j . 
An alternative, complementary understanding of the 
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FIG. 4: (Golor online) Paramagnetic phase diagram of degen- 
erate five-band model at half-filling. The dashed blue lines 
denote fixed ratios J/U between 1/10 and 1/2. The dots rep- 
resent the critical Coulomb energies of the Mott transitions 
along these lines. To the right (left) of the solid red curve, 
the system is insulating (metallic). The green box indicates 
the range of U and J values commonly used in DMFT studies 
of iron pnictides. For n = 6 the Mott phase occurs at much 
larger values oi U H, so that, approximately to the right of 
the solid red curve, the system is bad metallic. 



pseudogap can be achieved by realizing that the system 
at n = 6 occupancy exists in proximity to the n = 5 
Mott insulator. As shown in Fig. 4, in the limit of half- 
filling, the degenerate five-band model exhibits a Mott 
transition at a critical Coulomb energy that depends very 
sensitively on J, but is rather insensitive to the value 
oi U. In particular, Uc diminishes rapidly when J in- 
creases: f/c ~ 3 eV for J = 0.25 eV and C/c « 1 cV 
for J = 0.5 eV. (For J = C//4 and J = U/6, we find 
Uc w 1.5 eV and Uc « 2.0 eV, respectively, in agreement 
with Ref. 0-) In contrast, because of the high orbital 
degeneracy, J = yields a Mott gap only at very large 
Coulomb energy {Uc ~ 15 eV). Thus, for U and J val- 
ues usually employed in DMFT studies of iron pnictides 
(see green box), the half-filled system is a Mott insula- 
tor composed of localized S = 5/2 spins. Upon electron 
doping, some of the localized states become itinerant, so 
that the density of states exhibits a narrow quasiparticle 
peak below Ep, and a concomitant minimum or pseudo- 
gap above Ep, in agreement with the spectra shown in 
Fig. 2. DMFT calculations as a function of doping (not 
shown here) reveal that, away from half-filling, the upper 
Hubbard band quickly gets weaker and merges with the 
spectral feature above Ep, while the lower Hubbard peak 
remains rather stable. On the other hand, according to 
the phase diagram shown in Fig. 4, in the absence of ex- 
change interactions, the high orbital degeneracy yields a 
weakly correlated metal since U Uc- This also holds 
at finite electron doping, so that the spectral distribu- 
tion exhibits a broad quasi-particle peak, as indicated 
by the results for [/ = 3 eV, J = in Fig. 2(a). It is 
evident, therefore, that the collective mode in the self- 
energy shown in Fig. 3 and the associated pseudogap in 
the density of states (Fig. 2) obtained at realistic values 
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of J are a consequence of Hund coupling. 

As pointed out above, the complex geometry and band 
structure of pnictides gives rise to appreciable orbital 
dependence of the self-energy and interacting density 
of states, as has been found in various DMFT studies 

[l2-18[. In addition, screening processes differ greatly 
among pnictide compounds, leading to substantial varia- 
tion in Coulomb and exchange energies [l^l • Resonances 
in ImE(w) and the corresponding pseudogaps in the in- 
teracting density of states might appear therefore only in 
certain id orbitals. This seems to be the case, as indi- 
cated, for instance, by the kinks found in some Im E(iaj„) 
components of FeAsLaO Nonetheless, the robust- 

ness of the pseudogap in Fig. 2 over a wide range of U 
and J suggests that the characteristic peak-dip structure 
near Ep seen in quasiparticle spectra of pnictides ins- 
2ll . |26| might be partially due to Coulomb correlations. 
Optical spectra in the paramagnetic phase then should 
exhibit a depletion of spectral weight close to Ep, as was 
recently predicted theoretically |6| and observed experi- 
mentally in BaFe2Se2 0-01 ■ tti^ future, it would be 
very interesting to search more systematically for possi- 
ble links between pseudogaps and collective modes in the 
self-energy of pnictides, and to study their variation with 
doping and temperature. 

The picture discussed above has more general validity 
beyond the present five-band model. Thus, pseudogaps 
generated by self-energy resonances due to Hund cou- 
pling are also found in the analogous three-band system 
[301 . The main difference is that, because of the lower 
degeneracy, the boundary in Fig. 4 between metallic and 
insulating phases at half-filling is shifted to larger J (ex- 
cept at large U) [3l| and that at n = 3 ± 1 the system is 
relatively farther from the Mott phase than at n = 5 ± 1 . 

In summary, Hund coupling in ferropnictides plays an 
important role in the formation of a high-energy pseu- 
dogap. To disentangle this gap from density of states 
effects, the non-interacting part of the Hamiltonian is 
simplified in terms of five degenerate semi-elliptical sub- 
bands. The pseudogap can then be identified as a many- 
body feature that is stable over a wide range of Coulomb 
and exchange energies and that is related to a collective 
mode in the self-energy associated with spin fiuctuations. 
The nature of this mode can also be understood by view- 
ing the system as a doped n = 5 Mott insulator, where 
itinerant electrons coexist with localized spins. This pic- 
ture is consistent with the observation that at finite dop- 
ing the spectral distribution exhibits a narrow quasipar- 
ticle peak below Ep and a concomitant minimum above 
Ep, while most of the spectral weight resides in the lower 
Hubbard band. The crucial role of Hund exchange in this 
scenario is evident from the fact that, due to the high 
orbital degeneracy, at small values of J the pseudogap 
disappears and the system is a weakly correlated metal. 



Thus, the high-energy pseudogap is directly linked to the 
realistic magnitude of Hund's rule coupling. 

I like to thank Hiroshi Ishida for fruitful discussions. 
The DMFT calculations were carried out on the Jiilich 
Juropa machine. 
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